Tunneling Spectroscopy of Quasiparticle Bound States in a Spinful Josephson Junction 
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The spectrum of a segment of In As nanowire, confined between two superconducting leads, was 
measured as function of gate voltage and superconducting phase difference using a third normal- 
metal tunnel probe. Sub-gap resonances for odd electron occupancy — interpreted as bound states 
involving a confined electron and a quasiparticle from the superconducting leads, reminiscent of Yu- 
Shiba-Rusinov states — evolve into Kondo-related resonances at higher magnetic fields. An additional 
zero-bias peak of unknown origin is observed to coexist with the quasiparticle bound states. 



A quantum spin impurity in a superconductor modi- 
fies the order parameter locally via exchange interaction 
with quasiparticles, thereby creating Yu-Shiba-Rusinov 
sub-gap states [1-6 . For weak exchange interaction, a 
sub-gap state near the gap edge emerges from singlet 
correlations between the impurity and the quasiparticles. 
Increasing exchange interaction lowers the energy of the 
singlet state and increases a key physical parameter, the 
normal state Kondo temperature Tk- At ksTx ~ A 
(Kondo regime), where A is the superconducting gap, 
the energy gain from the singlet formation can exceed 
A, resulting in a level-crossing quantum phase transition 
(QPT) [6 -9 . The QPT changes the spin and the fermion 
parity of the superconductor-impurity ground state, and 
is marked by a peak in tunneling conductance at zero 
bias [10 . 

Yu-Shiba-Rusinov states can be studied in a controlled 
manner with a mesoscopic quantum-dot /superconductor 
Josephson junction [Fig. 1(a) and 1(b)] [11 . Such a de- 
vice provides a novel control knob that tunes the ex- 
change interaction — the superconducting phase differ- 
ence across the junction, ^. We label the states of 
the device \ndot,niead) , where Udot and niead represent 
the electron/quasiparticle occupation numbers in the 
dot/superconducting leads. When the spin orientations 
of the states are important, we replace the occupancy 
numbers with arrows. First, a spin 1/2 impurity is cre- 
ated by trapping a single electron in the lowest available 
orbital of the dot (assuming large level spacing) with a 
Coulomb barrier [Fig. 1(c)] p^il3j. This barrier is max- 
imal at the electron- hole (e-h) symmetry point, where 
the spinful state, |1,0), costs less than both the empty, 
|0, 0), and the doubly occupied, |2, 0), states by the charg- 
ing energy U. Charge fluctuation at energies below A 
is suppressed ioi U > A. Next, spin- flip scattering con- 
nects the degenerate states |t, \) and |^, t) via the virtual 
population of states |2,0) [Fig. 1(d)] or |0,0) [Fig. 1(e)]. 
These two scattering channels cause an effective (super-) 
exchange interaction between quasiparticles and the spin- 
ful dot. Finally, the exchange interaction is sensitive to 
phase due to the quantum interference of the spin-flip 
scattering amplitudes between the two leads. Compared 



to scattering via |2,0), scattering via |0,0) differs by a 
phase factor exp(— z0) because it is accompanied by a 
Cooper pair transfer [Fig. 1(e)]. At ^ = tt these two 
scattering channels interfere destructively, making the 
exchange coupling minimal at (/) = tt and maximal at 
(/) = 0. Consequently, both the singlet excited state, l^*), 
and the doublet ground state, \D)^ acquire a phase mod- 
ulation, albeit only in higher order processes for the lat- 
ter JJ^-TT. 

Previous experiments have investigated the ground 
state of spinful Josephson junctions [2'2H27]. Phase- 
biased junctions with weak coupling showed negative 
supercurrent [22| |23], consistent with theoretical pre- 
dictions of the weak phase- modulation of \D) p!4Hl6] . 
while for strong coupling, positive supercurrent was ob- 
served [24| ^ . The latter was interpreted in terms of a 
QPT associated with the interchange of states IS*) and 
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FIG. 1: (color online), (a), (b) Scanning electron micrographs of 
a lithographically identical device, (c) Lowest energy states of a 
single-orbital quantum dot at the electron-hole symmetry point for 
ksTx <C A. The states are labeled by their electron/quasiparticle 
occupation number in the format \ndot,niead) • Exchange inter- 
action dresses the states |1,0} and |1,1} as the doublet, |D}, and 
the singlet, states respectively. Transition from \D) to \S) 

produces a sub-gap resonance (SGR). (d), (e) Phase sensitive spin- 
flip processes coupling the |1,1} states \t,D and |^, t) via virtual 
occupation of (d) |2,0} and (e) |0, 0). 
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FIG. 2: (color online). Differential conductance as a function of 
tunnel-probe voltage, Vr, and back-gate voltage, Vbg- (a) Normal 
state data, B = 30 mT. (b) Superconducting state data, B ^ 
and (j) = 0. Coulomb diamonds in (a) and superconducting gap in 
(b) are highlighted with dotted lines. 



\D) at ksTK A p5H27]. Meanwhile, other experi- 
ments have performed tunnehng spectroscopy on spin- 
ful Josephson junctions without phase-control j28H3T] , or 
with phase-control but away from the Kondo regime [32] . 
This leaves the effect of phase on sub-gap states in the 
Kondo regime unaddressed. Tunneling spectroscopy in 
similar devices has also been used recently to examine 
signatures of Majorana end states [33H35j . 

In this Letter, we demonstrate both phase and gate 
control of sub-gap states in a Kondo-correlated Joseph- 
son junction {ksTx ^ A) [11]. We also report the first 
evidence of a singlet /doublet QPT induced by the su- 
perconducting phase difference. Our junction is an In As 
nanowire contacted by a pair of superconducting Al leads 
and a third normal metal tunnel probe. At magnetic 
fields above the critical field of Al, tunneling into the 
InAs quantum dot with odd electron occupancy showed 
Kondo resonances [12 with associated Kondo temper- 
atures, Tk ~ 1 K. Near zero field, tunneling into the 
nanowire revealed the superconducting gap of the Al 
leads, A 150 /ieV, and a pair of sub-gap resonances 
(SGR) symmetric about zero bias. Both gating the 
nanowire and varying the phase bias across the junction 
tune the SGRs. For certain parameters, the two SGRs 
cross at zero bias, which we interpret as a level-crossing 
QPT. However, no such crossing occured upon suppress- 
ing A to zero with an applied magnetic field. Instead, the 
SGRs evolve smoothly into Kondo resonances, and this 
transition is typically accompanied by the appearance of 
a separate zero-bias resonance of unknown origin. 

Our device consisted of a 10 fim long, 100 nm diameter 
InAs nanowire contacted by two ends of a superconduct- 
ing loop (5/100 nm Ti/Al). A third normal metal tunnel 
probe (5/100 nm Ti/Au) contacted the nanowire between 
the two superconducting leads [Figs. 1(a) and 1(b)]. The 
segment of nanowire between the superconducting con- 



tacts was 0.5 /im long, with a loop area of ~ 25 /im^. 
For this loop area, the flux period, /i/2e, corresponds 
to a perpendicular magnetic field period of 72 fiT. The 
nanowires were grown by molecular beam epitaxy on 
(lll)B-oriented InAs substrates using Au nanoparticles 
as catalyst, resulting in a wurtzite structure with a [111] 
growth direction. They were then liberated from the 
InAs growth substrate by sonication in methanol, and 
deposited on a degenerately doped Si substrate with a 
100 nm thermal oxide. Highly transparent barriers be- 
tween Al and InAs were made using an ammonium poly- 
sulfide etch immediately before deposition [36]. Lower 
transparency of the normal contact was achieved by etch- 
ing for less time. The device was measured in a dilution 
refrigerator with a base temperature 20 mK, through sev- 
eral stages of low-pass filtering and thermalization. 

When superconductivity in the entire device was sup- 
pressed by an applied magnetic field diamond patterns 
characteristic of weak Coulomb blockade (CB) were ob- 
served in transport between the loop and the normal lead 
[Fig. 2(a)]. Consecutive diamonds alternate in size, in- 
dicating that the orbital level spacing, ^, is comparable 
to the charging energy, U o:^ 200 /ieV. The smaller (odd 
occupancy) diamonds contain gate-independent zero-bias 
ridges that split at higher magnetic fields (Sup. 1), typ- 
ical of the Kondo effect [121 EZ] • From the temperature 
dependence of the zero-bias ridges, we estimate Tk to 
be in the range of 0.5-1 K (Sup. 3). Poor visibility of 
the odd diamonds suggests strong coupling to the super- 
conducting leads (r^ > /7), and the amplitudes of the 
Kondo ridges indicate an asymmetry between supercon- 
ducting and normal contacts [38\ The estimated asym- 
metry, Ftv ^ F^/IO, allows qualitative treatment of the 
Au lead as a weak tunnel probe. 

In the superconducting state {B ~ 0), gap-related 
features were observed at tunnel-probe voltages, Vt — 
±150 /iF iA/e, consistent with the gap of Al. SGRs 
symmetric about zero bias were also observed [Fig. 2(b)]. 
Comparison of Figs. 2(a) and 2(b) shows that the posi- 
tioning (in backgate voltage, Vbg) of SGRs in the super- 
conducting state coincide with CB and Kondo features 
in the normal state. In particular, the SGRs and their 
symmetric partners converge towards each other, to the 
extent that sometimes they overlap, in an odd CB valley. 
In contrast, they are pushed towards the gap edge in the 
even CB valleys. 

Based on their qualitative dependence on Vbg and 0, 
three categories of SGRs in the case of a spinful dot were 
identified: (i) For small charging energy, U < (A,r5'), 
SGRs do not cross the zero-bias axis for any Vbg or (j) 
[Figs. 3(a), 3(d), and 3(g)]. The SGR energy is maximal 
at ^ = and minimal at = tt [Fig. 3(d) and 3(g)] — this 
is the conventional phase dependence of non-interacting 
Josephson junctions [39]. (ii) For large charging energy, 
/7 > A, (Sup. 2) SGRs become inverted, crossing zero 
bias twice as a function of Vbg [Fig- 3(c)]. Between zero- 
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FIG. 3: (color online). Three sub-gap resonances (SGRs) arranged 
in columns of increasing Coulomb energy, (a)-(c) Back-gate depen- 
dence of the SGRs at = 0. The lower rows show their correspond- 
ing phase dependence off (d)-(f) and on (g)-(i) the electron-hole 
symmetry point, (d)-(g) Conventional phase dependence, (h) Hy- 
brid phase dependence, (i) 7r-shifted phase dependence 



bias crossings, the phase dependence of SGR energies is 
the opposite of the conventional behavior, that is, mini- 
mal at = and maximal at ^ = tt [Fig. 3(i)]. We call 
this a TT-shifted phase dependence. Outside the intersec- 
tions in Vbg^ the phase dependence of SGR energy is con- 
ventional [Fig. 3(f)]. (iii) For moderate charging energy 
U ^ A [Figs. 3(b), 3(e), and 3(h)], SGRs do not intersect 
for any Vbg at = [Fig. 3(b)]. Phase dependence away 
from the e-h symmetry point is conventional [Fig. 3(e)], 
but close to the symmetry point, the pair of SGRs inter- 
sect twice per phase period of 27r [Fig. 3(h)]. Grossings 
occur at (j) = tt ± ^0/2, where Scj) < tt is the phase- 
difference between the two closest crossings [Fig. 3(h)]. 
With this type of SGR, the phase dependence depends 
on the phase value itself: it is conventional for (p ^ 0, 
and TT-shifted for cj) ^ tt. 

Focusing on 7r-shifted SGRs, we examine in Fig. 4 the 
magnetic field evolution of three SGRs measured at their 
e-h symmetry points. The first SGR [Figs. 4(a)-4(c)] is 
the same as in Fig. 3(c). Selecting = from the full 
data set (Sup. 9), the well separated and resolved SGRs 
gradually approach zero-bias and merge into a Kondo 
resonance in the normal state [Fig. 4(b)]. Tempera- 
ture dependence of the normal-state Kondo peak gives 
Tk — IK [12] (Sup. 3). Taking ^ ~ 13 from normal- 
state GB data (Sup. 4), the splitting of the Kondo peak 
at ^ 140 mT is consistent with this value of Tk [40] 



[Fig. 4(b)]. In the other two cases (middle and lower rows 
of Fig. 4), Kondo peaks split at lower fields of 5 ~ 50 mT 
[Fig. 4(e)] and B < 20 mT [Fig. 4(h)], suggesting lower 
Kondo temperatures. 

In the second case [Figs. 4(d)-4(f)], SGRs overlap 
at zero-bias for (/) = 0, but are separated for </) = tt 
[Fig. 4(d)]. The overlapping SGRs at zero field evolve 
continuously into a Kondo resonance as the field is in- 
creased into the normal-state regime [Fig. 4(e)]. Phase 
dependent oscillations of the SGR vanish abruptly at a 
critical value of field. Be = 19.5 mT [Fig. 4(f)]. The same 
critical field is observed in Fig. 4(c), and also in higher 
density regimes of the device (Sup. 5). 

The last case in Fig. 4 has no phase-dependence 
[Fig. 4(g)], presumably because of poor coupling to one 
of the superconducting contacts. However, its Vbg de- 
pendence allows us to establish that this SGR is indeed a 
TT-shifted type (Sup. 8). Here, in contrast to the first two 
cases, the pair of SGRs evolve continuously and directly 
into a split Kondo peak without ever merging or crossing 
at zero bias [Figs. 4(h) and 4(i)]. 

Glose inspection of Fig. 4 reveals an unexpected and 
intriguing feature: a narrow needle-like resonance pinned 
at zero bias. In Fig. 4(b), this "needle" is absent at 5 = 
but appears for B > 10 mT while the leads are still su- 
perconducting. In Fig. 4(d) the needle is hidden by the 
SGRs at = 0, yet it is clearly visible at = tt. In this 
case, the needle exists at 5 = 0, and merges into the 
normal-state Kondo resonance at higher field [Sup. 9(f)]. 
In Fig. 4(h), the needle appears at 5 > 10 mT, simi- 
lar to the case in Fig. 4(b), despite a large difference in 
Kondo temperatures. In fact, the strength of the needle 
appears uncorrelated with Tk of the normal-state Kondo 
peak (Sup. 10). The needle is also distinct from the nor- 
mal state Kondo resonance as seen in Figs. 4(h) and 4(i), 
where three separate peaks can be identified: The two 
peaks fianking the central needle appear to emerge from 
the SGR at the low-field end and evolve continuously 
into the split Kondo peaks at the the high-field end. We 
find that the needle only appears between the two Vbg 
intersection points of 7r-shifted SGRs, which in turn cor- 
responds to an odd Goulomb diamond (Sup. 6). Finally, 
the needle appears brighter at = 0, when the separa- 
tion between the two SGR is the smallest [Figs. 4(c) and 
4(d)] (Sup. 6). 

We now compare theoretical expectations for 
SGRs [17 to experimental observations. At the 
e-h symmetry point of a spinful quantum dot with sup- 
pressed charge fiuctuations, the phase-tunable exchange 
interaction detaches a singlet state \S) down from the 
gap edge [Fig. 1(c)]. Since quantum interference weakens 
the exchange interaction at ^ = tt [Figs. 1(d) and 1(e)], 
a TT-shifted SGR is indeed expected (phase modulation 
of the energy of \D)^ being a higher-order effect, is much 
weaker than that of l^*)) [181421] . This is consistent with 
out experiment, as seen, for example, in Fig. 3(i). Strong 
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FIG. 4: (color online). Arranged in the order of decreasing Tx, each row shows the evolution of a SGR at the electron-hole symmetry 
point as a function of phase and magnetic field. The left column shows phase dependence at 5 ~ 0, the center column shows magnetic 
field dependence at = 0, and the right column shows the magnetic field and phase dependence around B = 18 mT. To obtain the phase 
constant panels (b) and (e), we select (j) = data points from the full data set. The oscillations of the SGRs disappear abruptly at B = 
19.5 mT (dotted lines) in both (c) and (f). Inset in (b) is a closeup of the region outlined with dotted lines. A third resonance, pinned at 
zero bias, is clearly visible in the high contrast color scale. 



coupling to the leads, reflected in the large T^, should 
further result in a SGR that is well separated from 
the gap edge at = ^ [9]. Detuning Vbg towards a 
neighboring even diamond increases charge fluctuations 
and mixes either |0,0) or |2,0) into l^*), thereby lowering 
its energy. Consequently, one expects a level-crossing 
QPT to a singlet ground state as Vbg approaches an 
even diamond, in agreement with the zero-bias crossings 
in Fig. 2(b) and Fig. 3(c). This QPT is predominantly 
governed by the enhanced charge flutuations away 
from the e-h symmetry spot. Finally, the observed 
conventional phase dependence in the even state of the 
dot [Fig. 3(f)] is also expected, because a spinless dot 
acts purely as a potential scatterer, behaving as an 
effectively non-interacting junction [41 . 

A more interesting QPT occurs in Fig. 3(h) as a func- 
tion of phase-bias. It corresponds to a situation, where 
the energy gain from the quasiparticle-dot singlet forma- 
tion makes this state the ground state at ^ = 0, but not 
at ^ = TT. This behavior is known in the theory literature 
as -junction or tt -junction [l7l[42], and, to our knowl- 
edge, has not been reported in previous experiments. 

Reducing A sufficiently below ksTx should result in 
a level-crossing QPT that is driven entirely by spin fluc- 
tuations [6]. Experimentally, we would see a zero-bias 
crossing of the SGRs at 5 < 5c as 5 is increased to 
suppress A. However, this theoretical expectation is not 
seen in our device as exemplified in Figs. 4(b) and 4(h), 
perhaps obscured by our current experimental resolution 
or by the needle feature. The needle may be related to 



similar features observed in recent experiments ^Ql [31] . 
An unlikely soft gap in Al may explain such a resonance 
in terms of conventional Kondo screening. We note, how- 
ever, that the needle itself does not split with increasing 
5, as one might expect from a conventional Kondo effect. 
More intriguingly, the needle appears much stronger at 
(/) = than at (/) = TT, suggesting possible phase depen- 
dence and a link to the sub-gap states (Sup. 11). While 
the observed behaviors of sub-gap states agree at 5 ~ 
with existing theory on Yu-Shiba-Rusinov states, further 
theory and experiment are needed to understand the ori- 
gin of the needle and the magnetic field dependence of 
the sub-gap states |43l |44] . 
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